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ABSTRACT 


A unified approach is applied to the derivation of a number of formulas 
for the probability of signal detection and the probability of false alarm. 

The context is incoherent integration, with fluctuating signal-to-noise ratios 
and/or fluctuating thresholds. The standard results are obtained and extended 
to more general fluctuation models. A fundamental duality is established 
between fluctuating signals and fluctuating thresholds and used to simplify 
the derivations. Also included is an expression of the cumulative 
F-distribution as a finite sum of llarcum Q-functions. 
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I. INTRODUCTION 


In a wide variety of signal detection and binary decision problems, the 
evaluation of performance reduces to the computation of a number, P, defined 
as the probability that one random variable, u., exceeds another, called y. 

In symbols, 

P = Frob {u > y} , 


and u represents a processor output which is being tested for signal preserve, 
while y is a threshold, which may be a constant. In all cases considered 
here, u and y are independant, and u is the result of incoherent summation of 
N complex-samples each containing signal and noise components. These 
quantities, z^, represent I and Q output samples of some coherent processor 
and they are modelled as sums of signal and noise terms: 

z = s + w . 
n n n 

The noise components are Gaussian, independent and with zero-mean, all sharing 
the same variance: 

E|w I 2 = 2o 2 , 

n 

and all having the circular property 

Ew 2 = 0. 
n 


Various models will be used for the signal components, corresponding to (and 
extending) most of the "fluctuating target" models of the radar literature. 
The random variable u is the normalized sum 


l N 

u = —s- f Iz | 
2„ 2 nil n 


( 1 - 1 ) 


1 











hence its probability distribution function (pdf), conditioned on the signal 
sequence, {s n ) , is 


f N (u,a) = e U a (u/a) 1^ 


( 1 - 2 ) 


In this formula. 


a 



N 


l 

n~l 



(1-3) 


and I n is the Bessel function of imaginary argument. Since the signal 
components enter this conditioned pdf only through the sum a, the "random 
signal" models will be represented by postulated pdf’s for a. 

In a similar way, "fluctuating thresholds," such as those arising in CFAR 
problems, are represented by putting y=cx, where c is a scalar multiplier, 
and x is a random variable. Various pdf's are postulated for x, including the 
pdf f^(x,d), with fixed or fluctuating parameter, d. Ordinary "linear CFAR” 
is represented by this model with d=0. 

'2 3) 

When a and y are constant, P is the standard Marcum Q-function^ ’ ' which 
must be computed from an infinite series, used together with a bound on 

( 4 ) 

truncation error. An efficient algorithm has been described by Shnidman, 
which involves recursive computation of the terms, along with several 
refinements which improve computational efficiency. The standard series, 
(Fehlner's formula, see Section 6) was used by Mitchell and Walkeras the 
basis of a systematic derivation of expressions for some of the cases 
presented here. Because they start with an infinite series for the Q-function 
and obtain their results by averaging over various random models for a and/or 
y, all the expressions obtained are in infinite series form. Although all 
these series solutions lend themselves readily to recursive computation, it 
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happens chat in most of the cases (except the original one, where both a and y 
are constant), finite-sum expressions exist and many of these can be found in 
the literature.^ 

In this study, we start with an integral representation for the 
Q-function, instead of an infinite series. Then, in analogy to the method of 
Mitchell and Walker, we derive the desired results for random a and/or y by 
averaging this expression under the sign of integration. Then, by simple 
changes of the variable of integration, the integral representation is itself 
used to obtain the desired finite-sum expressions. In all but two cases, the 
terms of these sums are elementary functions, easily computed recursively. In 
the two exceptions, the terms of the sums are themselves Q-functions. One of 
these, actually representing the cumulative F-distribution, is the source of a 
series of generalizations of the standard fluctuation models. 

When a=0 in formula (1-2), the resulting pdf will be called g^(u): 


M JL 

g N (u) « f N (u,0) e 


0-4) 


This function, the chi-squared with 2N degrees of freedom, is used to 
model signal and threshold fluctuations for the standard formulas derived in 
Section 4. Only integral values of N are used, since non-integral values (as 
in the Weinstock cases^^) do not lead to detection formulas of finite form. 
Fluctuating signals are modelled by postulating that the pdf of a is 


V a > * I %<t> 


(1-5) 


where M and b are fluctuation parameters. In this model, the mean value of a 


a = Mb, 

and according to (1-3), 


( 1 - 6 ) 














1 N 7 

I = -V l < s i 

2c 2 n-1 " 


(1-7) 


We denote the average signal-to*-noise ratio for a single sample by “SNR,” so 
that 


SNR = | 


for a fixed signal, and 

SNR = | = | b (1-8) 

for a signal fluctuating according to (1-5). The standard “Swerling cases” 
correspond to M-values, as follows:' 2 ^ 


Swerling Case M Value 

1 1 

2 N 

3 2 

4 2N 


When (1-4) is used to describe a fluctuating threshold, we write 
y = cx, 

and use 

g L ( x) 

for the pdf of x. The pdf of y is then 


f (y) = — 

o c 


*L<t> 


(1-9) 


4 










This model represents linear CFAR, and the variable x is the result of 

incoherent summation of L complex samples containing only noise of the same 

variance as that of “he signal samples. The terms of this sum are normalized, 

2 

just as in (1-1), by the true noise variance, a , but these normalizing 
factors cancel when the probability that u exceeds y is computed. If the 
noise samples are known to have a different variance than the signal samples, 
the variance of the signal samples is used for normalization and the variance 
ratio is absorbed into the constant, c. Note that x is the sum of the noise 
samples, not the average, which affects only the significance of c. 

The fluctuation models are generalized in Section 5. First, the pdf, 

£ o (a > -H ( t- d) - w - io) 

is used to describe signal fluctuations. The expected value of a in this case 
is a = (M4d)b, hence 

SNR = | b(l + £) . (1-11) 

This case is of rather academic interest, and the dual problem with fixed 
signal and threshold y=cx, is more interesting. The pdf of x is taken to be 

f Q (x) = f L (x,d), 

which models linear CFAR, with L samples in the threshold, but with signal 
components (unwanted in the CFAR application) included in these samples. Only 
the sum of these signal powers affects the pdf, and the average 
signal-to-noise ratio, per threshold sample, will he d/L. 

In further generalizations, a is allowed to fluctuate, using the same 
models as in Section 4 and also (separately and together) the new threshold 
signal components are randomized using the pdf 





















( 1 - 12 ) 


f o< d > - K H ( b - 


for the parameter d. The average SNR, per threshold sample, in this last case 
is Kh/L. 

A key feature of the basic approach is the duality between fluctuating 

signals and fluctuating thresholds. In Section 3 a useful artifice is 

introduced which allows us to obtain the detection probability for a fixed 

signal and threshold fluctuating according to a pdf, f Q (y), from the solution 

of the dual problem, where y is constant and a fluctuates with the same pdf, 

f (a), 
o 

The results are collected without proof in Section 2, and the basic 
method is developed in Section 3. The standard formulas are obtained in 
Section 4, and some extensions of these results are given in Section 5. 

Finally, in Section 6, some alternative finite forms are obtained and 
for ore common case a derivation of the finite fora directly from the 
equivalent infinite series is given. 

n. COLLECTED RESULTS 

The results of this study are collected here without proofs, all of which 
will be found in the next two sections. In all cases, the assumptions and 
conditions are spelled out in detail, but detection and false alarm 
probabilities are called simply and P„ A , since the number and variety of 
parameters become too great to list as arguments in some of the cases. The 
following conventions are used for the most coimnon parameters: 


Parameter 

N 

L 

M 


Significance 

Number of "signal-plus-noise” samples 
integrated for detection. 

Humber of "noise-only" samples 
Integrated to establish a CFAR threshold. 

Signal fluctuation parameter, corresponding 
to chi-squared distribution with 2M degrees 
of freedom. 


6 






K 


Fluctuation parameter for signal components 
in threshold samples, again chi-squared, with 
2K degrees of freedom. 

a Signal parameter for fixed signals 

b Signal parameter for fluctuating signals 

c CFAR threshold multiplier 

d Parameter for fluctuating signal 

components in threshold samples 

Three functions enter these formulas. These are defined here and 
discussed elsewhere. 


(1) P N (y,a) = / f„(u,a) du , (2-1) 

y 

where f,. is given by (1.2). P„(y,a) is the probability of detection for 

i» a 

fixed, normalized total signal-to-noise ratio, a. It is essentially the 
Marcum O-function for incoherent integration of K samples. 


(2) s s (,) - P K (y,0) - e'* Y k 

m=0 


( 2 - 2 ) 


S,,(y) is the false alarm probability corresponding to ?„(y,a). It is computed 

a a 

recursively: 


s,„, (y) = s,.(y) + T„(y) 


1 
















N-l 


(3) R N (y»M) = l ( M "J; +n ) y n . 


n=0 


These functions are also computed recursively: 


R N+! (y,M > = V y ’ M) + V y * M) 


V N (y,M) - (1 + ~)y V N _ 1 (y,M) / 


R^y.M) = V q (y,M) - 1 


N>1 


(2-4) 


(2-5) 


and 


The probability density functions used for the flu -uation models are 


N-l 

%<*> ■ (fejr ^ 


N-l 

2 


f N (x,a) » e (-) 1 n _ 1 (2/ox) 


A. Non-Fluctuating Signal 

In these formulas the signal components of the random variable being 
tested for signal presence are constants. The average SNR per sample in all 
these cases is 

SNR = a/N. 


(1) Fixed Threshold, y 


8 
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(4) Fluctuating Threshold, y=cx, conditional pdf of x is 

9 









f (x|d) ■ f (x,d) ; pdf of d is 

O L. 


f 0 (d) " ¥ g K ( f ) * 


Condition: K>L 


p „ ! _ c N (l+h) N+L ~ K K y L ( K-L ( __h_ ) J 

(l+c+ch)^ -1 iio * 1+C+Ch ' 


* Vl <I*£h> 

n=0 


p . (l+h) L ~ K __ _ K y L ( K-L ( h * 

FA d^h)^- 1 tio * 1+c+ch 


V ,N+L-l+£v , . , .n 

x l ( ) (c+ch) 


B. Fluctuating Signal, pdf of a is 


V’ “ b 8 M^ b ^ 


Average SNR per sample: 


(1) Fixed Threshold, y, P„, given by (2-7) 


(2-12) 


(2-13) 


10 








M>N: 


P D - 


Cl+b) 


1 M " N U-M . 

< A > b S N+A ( l+b> 


M<N: 


P D “ 


(^) N - M [? 1 («:* 


'> ( - b) ‘ Vife) 


N-M A | 

£ <Vl‘> <ife> V y) | 


Special cases: 


M-N: 


P D " S N ( I+b ) 


M-2N: 


P D " 


ci+b) N io < ‘ ) bt Swt 


M-l (N>2): 

1+h 

• [i 


' S N-1 ( & )] 


M-2 (N>3): 


l+b, N " 2 " 1+b 


P D * S N-2<*> + <T> • 


|a - N b 2 ) u - s N _2 ( iTb )J 


(2-14) 


(2-15) 


(2-16) 


(2-17) 


(2-18) 


11 



... . .....pip*................... 


+ 


y 

l+b 


(2-19) 




M*2N: 


l (") b 1 . « 


N 


D (l+b+c) L 


. a 

■ V 0 “N+-£ v 1+b+c 


(2-23) 
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ti»l: 


(1+b) 


N+L-l 


N-l 


D b N - 1 (l-.b + c) L 


( He ) 1 *"' 1 £^0 ' 4 


" 1 ‘ i+crrc 


M=2: 


(14,b)N ^. L - 1 - jl - (N-2) ^ + (N+L-l) C 

L+l b 


D b N_2 (l+b+c) 


1+b 


w 4 - ii+<£+2 - N)i ^ 1 ^ 

) 1=0 ( 


(1+c) 


£+ 2 - 


(3) Fluctuating threshold, y=cx, pdf of x is 


f (x) = f (x,d). 
O 


P_, given by (2-11) 
FA 


M>N: 


D d+b+c ) 84 - 1 -' 1 A * 1+b+c 


x 


N+Z-l 

I 


n=0 


N+L+ £-1s,1+bs c , cd ^ 
■ L+n K c } n+l4+b+c ; 
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Special cases: 


M=N: 


„ _ c N ^l+b) 1 V ,N+L-l x/ l+b, n e , cd _ x 

P ~- IT ^* 1 { L+n K c ; n+ri+b+c ; 


(Ifb+c )^- 1 


n=0 


11=2K: 


p . 

D ' 


*■ V C l+b+c' 


£=0 


N+£~l TIL 11 J 

r / N+L+?-1 n , l+b N „ , cd x 

* ' L4r. )( —> WltS*’ 

n"0 


li*l: 


P = (l+b)^" 1 1+c+d 

D “ b^dd-b^) 1 ^ 


N-l H-l M , T , 

c i x 


„ . X N+L-1 ^ ' L+£ 

(1+c) £=0 


cd 


, ,l x £ o , cd x 1+c+d ,l+b+c x ,, r bed 

x * ( c ) s £+i ( t+c } ■ 6 ( -bT~ ) S £+1 1 a+eTcr+b+c) 


(4) Fluctuating threshold, y=cx, conditional pdf of x is 

f (xid) ■ f T (x.d) , pdf of d is 
O Li 

V d > ■ ¥ 8 k<¥> 

Any H, K>L: 


(2-27) 


(2-28) 


(2-29) 


14 















, cVjj. ) atL - K (14c4.h) H - |< - 1 - +1 x 

P = 1-- 

" (1+b+c+ch) 


K ; L ,K-L V h * L r £_1 .IH-L+t-iw ^ h) n R ( b , M) (2-30) 
I < « )( 1+r+ 7h ) L ( N+n Hc+ch; n+l k 1+b+c+ch 


* £-0 ' * /Vl+C+Ch n=0 


Pp A given by (2-13) 


Any K, M>N: 


P D' 




(1+b+c+ch) 


K 


M ; N ,M-H W be * N+ r _1 ^L+A-lwi+k^ r (- - , K) (2-31) 

x I ^ e ><T+b+T> ' ' L+n ^ c ' n+1'‘1+b+c+ch 

£=0 n=0 


N-l n K-N-L+ 1 N-l 
c (1+c) 


'FA 


(1-fc+ch) 

Special Cases: 

K=L, any M: 


K 


r x ,N+L-1. ,1 » „ / ch R v (2-32) 

X < L+n } ( c> R n+l l+c+cir’ K) U ^ 


n=0 


. AihwW)"^ 

P = 1 - --- w 

D (1+b+c+ch) 


L-l 


V* .N+L-lx, ^.xii „ , b 

I ( w+n R n +1 1+b+c+ch 


M) 


n=0 


(2-33) 


15 
















(2-34) 


P 


FA 


1 


.N+L-l 

1 . 1-1 c-i-cn; 


N—1 

I ( 

u— 0 


N+L-l 

n 


)(c+ch) n 


I1=N, any K: 


P 


D 


N-l 


(1+b) (1+b+c) 


K-N-L+l 


(1+b+c+ch)^ 


N-l 

l 

n=0 


. N+L-l v 1+b' 
1 L+n n c ; 


n 


R ^ ( 


ch 


h+1 1+b+c+ch 


K) 

(2-35) 


P FA given b y (2-32) 


C. Fluctuating signal, pdf of a is 


f (a) = i f M (~ 
o o a b 


d) . 


Average SNR per sample: SNR = ^ (!l+d) 

(1) Fixed threshold, y. Pp^ given by (2-7) 
ii >N: 


P D = 


(1+b) 


M-N 


ti-N 

l 

£=0 


(li-N) 
v l ' 


b P. 


( 


y 


bd 


N+£ v l+b ’ 1+b 


) 


(2-36) 


It is often required to evaluate detection performance when the false 
alarm probability is fixed, which involves the inversion of the formula for 
Pp A in order to determine the threshold. The Newton-Raphson iterative 
method(^) is useful for this purpose, and we discuss its application to the 
two most common false alarm formulas, (2-7) and (2-9). The other cases listed 
above involve CFAR problems with signal components in the thresho d samples, 
but these components are usually not anticipated, and the threshoi. will have 
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been found from (2-9). Formulas like (2-11) and (2-13) then show the effect 
of the signal components on Pp A after the threshold has been fixed. 

Beginning with (2-7), we want to solve the equation, 

$(y) 2 S N (y) - p pA - 0 , 

for y, where P pA is the assigned probability of false alarm. The 
Newton-Raphson iteration is 


y 


n+1 


0<v n ), 


where 


G(v) 


$(y) 

$'(y) 




and y c is a suitable initial value. A good plan is to start the iteration 
at a point where the derivative, $’(y), is large, to avoid wild oscillation of 
the sequence, y^. This can be assured by choosing v q to be the solution of 

♦-(y Q > = o. 


which will be unique in the applications made here. 
From (2-2) we have 


$’(y) “ S N (y) = -e 


-y 


N—1 

m 

N-2 

m 

V 

Ve'" 

m! 

y 

y_. 

a 

ii f 

o 

m=0 

m! 
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In C2—3) we nave written, in effect, 


N-l 

s„(y) - I 

n=0 


where 




» 


hence 


S N (y) = "Vl (y) * 


As Sy(y) is computed recursively, from (2-3), for use in the iteration, the 
final term, T^_^(y), can be recovered for the evaluation of $'(y). Then 


G(y) 


y + 


s N ( y ) - p fa 

W y) 


> 


and iteration is continued until the correction term, G(y)-y, reaches some 
preassigned small value. We also find 


r(y) 



\1 


N-l 

[y -(N-l)y 



> 













and so 


y = N—l 
o 


is the desired starting value. 


A very similar procedure can be used for (2-9) with 


1 r , N+L-l s _.£ „ 

4>(y) o^nN+L- 1 - * £ 5 * " ‘FA* 

(1+y) £=0 


(we have retained the symbol y, instead of the c of (2—9)). 
This time. 


\ N+L-l r ,m-l, . £ 

» ~ ~ ■■ .IM. .1 ' l ) y 

(1+y) £=0 


N—1 

, 1 V (N+L-l)? £-1 

( l + y )^ 1 ’ 1 £ 1 (N+L-£~x) !(£-!)! > 


-^lY C*^ 1 ) y' - Y (■£*><,« ♦ yS 

(l+y)‘ L < £=0 £ £=1 1 1 


Using the reduction formula 


, N+L-l, = -NbL-2. ,N+L-2. 

^ £ ’ ' £ } K £-1 1 ' 


the quantity in curly brackets becomes 



... ii[ i|(iip i« win ii‘ iiimii +* iiiilll|lll||lll|||i|||iilHIIIIIIIIIII”ill |l ^ll l ii 




















N-l 


UiT O 


N-l 




C * r. ,HTL-4, X 

1 ( c ) y - > ^ , ) y 

£=0 X £=1 * 1 


»»*] *» 
N-i } y 


Therefore 


♦*(y) = - 


N+L—1 ,N+L—2» N—1 

■< ) y 


(l+y) 


N+L 


L . 1 .N+L—1- N-l 

1+y Cl - K ')^" 1 ! ' _1 7 ' 


and if we define 


n fy) = _1_ "y 1 „* 

V 4-0 ‘ J 


N-l 


= I Q»(y)» 
£=0 


then 


q £ (y> 


1 ,N+L-1» £ 

< * > y 


CI -^- 1 * 


and 


*'(y) 


-1% Vl^- 


Thus the iterative algorithm is based upon 
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G(y> = 



w N (y) - p 


FA 


5 


and , (y) is recovered 
it—i 


fron the iterative coaputation of W (y): 


w (y) 

XI 


Vl<*> + W y) 


,(y) = - l)y Q £ _ 2 (y) 


^(y) - Q 0 (y) 


(i-*y> 


N+L-i 


Finally, we compute 


/U+L—1, (K-l)y K ~ 2 (H 7 )-(lfrt.)y" 1 
N “ 1 


and therefore 


v 

'O 


L+i 


provides a suitable initial value, 
in, BASIC FOHHULATIOS 

Tna whole analysis here is based on Schlafli’s integral 

( -Q\ 

't for the Bessel functions. We use it in the form 
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I n (x) 


n i 
v 2 ; 2 tt1 



2 

x 

4t dt 


t 


n+1 


( 3 - 1 ) 


The contour of integration is a small circle, enclosing the origin in a 
positive sense. With this representation, the pdf of u (for fixed a). 


N-l 


f N (u,a) = e u a (J) 2 1^(2®, 


( 3 - 2 ) 


can be written 


t + 


au 


c , s -u-a N-l 1 r t dt 

V“' a) = e “ sr k |. e e 7 


Replacing t by ut in the integral, we obtain 


f N (u,a) = e 


-u-a 1 
2iri 


ut + — 

t dt 


|tl=e 


N * 


( 3 - 3 ) 


This representation is more convenient than the Fourier transform 
expression., which can be obtained from it by the. change of variable t=l-iX. 

We carry out this transformation, since we need the characteristic function of 
f^ later or.. After substitution, we have 


f N (u,a) = e 


1 

2ni 


-iuX + 


/ e 
C 


a 

1-iX 


(-idX) 

(l"iX) N 
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where C is a small circle enclosing the point X * -i. Tills contour can be 
expanded, without changing the integral, until it runs along the real X-axis 
from +A to -A, and then back along an arc of radius A, in the lower half- 
plane, until it closes at X * +A again. Because of the factor exp(-iuX) in 
the integrand, the integral along the arc vanishes in the limit A+«, for any 
N=l,2,.... We then reverse the path along the real axis to obtain 


f N (u,a) = e 


00 


/ 


-iuX 

e 


a 

1-iX 

e 

(l-iX) N 


dX 

2tt 


The inverse transform gives the characteristic function 


« . „„ 1-iX 

/ e f„(u,a)du = e~ n -- = fc/X.a) • 

0 N (1-iXr 


(3-/.) 


Recalling the definition of a in equation (1-3), we can write <j> N (X,a) as a 
product: 


¥ x) 


N 

7T 

n=l 


-a 


n 


_n 

1-iX 

e 

1-iX 


2 2 

where a n =• Is^l /2a , which corresponds to the definition of the random 
variable, u, as a sum. Each factor here is the characteristic function of a 
non-central chi-squared random variable, with two degrees of freedom, and u is 
non-central chi-squared with 2N degrees of freedom. 

From our basic integral representation for f^(u,a) we obtain the "no 
signal" special case, 
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This set of pdf’s will be used to model fluctuating thresholds and signal 
parameters later on, and we introduce the notation 


S N < U > = V u >°) • 

Evaluating the contour integral for f^(u,0), we get 

/ X _ 1 N-l -u 

8 N (u) “ (N-l)! u e * 


( 3 - 5 ) 


which is, of course, the chi-squared pdf, for 2H degrees of freedom. 
Obviously, Eu«N, if g^(u) is the pdf of u. 

The "probability of detection" corresponding to the pdf f N (u,a) is 
denoted ?^(y,a) and it is defined by the equation 


P N (y»a) = / f N (u,a)du . 

y 


( 3 - 6 ) 


Our notation differs from the standard one (4) in that we retain the same order 
of the variables as in the pdf, and the parameter a is the "total SNR," as 
defined Ln Section 1. The advantages of this choice seem worth the cost of 
departing from common usage. In terms of P„, Marcum's: Q-function is 













To obtain an integral representation for P^, we substitute the 
representation (3-3) in definition (3-6) and reverse the order of 
integration: 


P N (y,a> - e' a 


( t 

It - 8 


t -(l-t)u dt 

J e du 17 


Since Re{l-t} > 0 on the t-contour, the u-integral is uniformly convergent, 
thus validating the interchange, and we obtain 


? N (y,a) = e 


-y-a 1 


yt + t dt 


111 =*e 


t N (l-t) 


(3-7) 


This is the fundamental representation for P^, used throughout this study. 

A second representation is obtained by making the change of variable, t+l/t, 
in (3-7): 


p N (y»a) - e 


-y-a 1 
“ 6 2iri 


I 11 =1 / G 


* t+ £ t^dt 

5 t-1 


The contour now is a large circle, including the simple pole at t=l, as well 
as the essential singularity at the origin. The residue of the integral at 
the simple pole is exp(a+y), and hence we can write 
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„ , -y-a 1 , at + t t N_1 dt 

P^Cy.a) - 1 - e Til , e "Ft 

I t|=e 


(3-8) 


where the contour is again a small circle enclosing the origin. 

Both representations are useful and we note that the integral in (3-8) is 
just like the one in (3-7), except that the variables y and a are 
interchanged, and the factor t appears with a non-negative exponent. In other 
words, we have the evaluations 


-y-a 1 


yt + t t" n dt 


? n (y,a) ; n>0 


1-P ln (a,y) ; n<0 


It proves to be exceedingly useful to define P N (y,a) for all integral N by 
means of representation (3-7), with the understanding that 


P_»(y»a) si- P (a,y). 


(3-9) 


This artifice allows us to exploit the duality between threshold and 
signal parameter In the following way. Suppose the signal parameter is 
randomized, according to some pdf, f Q (a)* With fixed threshold, y, the 
probability of detection will be 


F M (y) = I P M (y,a)f (a)da. 


(3-10) 
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Now suppose the right side is evaluated for positive and negative N, using 
(3-7) as the definition of P^(y,a). Then the dual problem, in which a is 
fixed and y is random (with the same pdf) has detection probability 


00 CO 

/ P M (y,a)f (y)dy = / U - P „(a,y)] f (y) dy 
0 N ° o 1 N 0 


= 1 - Fj_ N (a) • (3-11) 

This technique saves the needless duplication of calculations, which would 
otherwise use (3-7) for one set of problems, and (3-8), in a parallel way, 
with the dual problems. 

The source of this duality is, of course, the basic pdf. The fact |hat 
l n (x) - I_ n (x) can be proved from (3-1) by making the substitution t 
in the integral and hence from (3-2), 


W u ’ a) = e 


-u-a 



-a-u 

e 



N 

2 F^C^au) 


Thus we could define 


f l-N (a,u) = f W4-i (u ’ a) * 


m' 


(3-12) 
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and this definition is consistent with (3-6), extended to negative values of 
N. Integration of both sides of (3-12) over a yields 


/ f 1 _ N ( a » u > da 3 P i-N (b ’ u) = 1 " V' u » b) 

b 


° / f N+ i (u,a)da ’ 

b 

or 

00 

P N (y,a) = 1 - / f N+1 (y,a')da', (3-13) 

a 


an interesting result which also follows directly from (3-3). 

The probability of false alarm for fixed threshold will be denoted 

S N (y) = V y ’ 0) * 


From (3-7) we have the integral representation 


S N (y) 3 ^ 2¥1 


[ e yt —gi¬ 
lt I =e t N (l- 


t) 


(3-14) 


which will often be used to ivaluate contour integrals in this analysis. For 
N>1, we have 


00 
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after expanding (1-t) \ and then 


N-l n 

S M (y) = e ^ £ ^~r • 

N J L n m! 

m=0 


(3-15) 


The series terminates, since there is no pole at t=0 when m>N. Of course. 


s N - (y ) " I S N (u)du , 

y 


which again yields (3-15) after substitution of (3-5) and repeated integration 
by parts. 

From (3-14) we have 


S (x) ® 0 , n<0 . 

n 


(3-16) 


This is completely consistent with (3-9), since 


S_ N (x) 3 P _ N < X *°) 3 1 ~ P n+ 1^°’ x ^’ 


and 


V°» x > 


3 I 


f N (u,x)du = 1, 


by normalization. 

We proceed now with the general problem of detection with a fixed 
threshold and random signal parameter. As above we let f Q (a) be the pdf of 
a, and we define the characteristic function 
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(3-17) 


j> (X) = / e^ a f (a) da 

0 0 0 


We evaluate the detection probability by substituting the integral 
representation, (3-7) into Eq. (3-10): 


/ £ o <a) 


-y-a 1 
s 2iri 


I t{=e 


yt + f dt 


t N (l-t) 


da 


We want to reverse the order of integration, which will involve the evaluation 
of 


-(1 - 


I )a 


f (a) da . 
o 


For this integral to converge uniformly, it is necessary that Re(l- t) > 0, 
which will be satisfied if we can arrange to have !11 >1. This is accomplished 
by expanding the original contour. It I = e, until it becomes a circle with 
radius larger than unity, and compensating for the for the effect of the pole 
(at t=l) now included. In other words we write 


-y-a 1 
2ui 


|t|=e 


yt + - 


dt 


t (1-t) 
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= 1 + e 


-y-a 


1 

2iri 


ltl>l 


yt + f 


dt 


t N (l-t) 


and Chen reverse the order of integration. The result is 


F N (y) = 1 + e 


2ni 




[ t| >i 


00 -(1- ■—)a 
/ e f (a)da 

is o 


dt 


N., , 

t <l-t) 


The integral over a yields $ [i(l- ■£•)}, with uniform convergence to justify 
the interchange. 

Now 


U D (X)I < 1 

whenever Im(A) > 0, and this upper half-plane maps into the exterior (and 
boundary) of the circle. 


It 



1 

2 » 


in the t-plane. It will be true in the cases considered here that 4> 0 (A) is 
analytic when Im A > - e, for a suitable positive e, so that <j> o (i(l - ^-)] 
will remain analytic for a small but finite distance inside the circle 

It --^1 “ j. This will allow the t-contour to be shrunk just inside the unit 

circle. Since <f> Q (0) = 1, the pole at t=l is simple, and hence 
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dt 


F (y) = 1 + e' y f e yt 0> [i(l - - 

ivl |tI>1 0 C t (1-t) 


or 


V*> 


-y 1 r yt , rw , l N , dt 

* e ■=—r j e $ [i(l - -)] —- 

2tti if, o t N,, . 

ItI=l-e t (1-t) 


(3-18) 


This is the basic expression which will be used throughout this report. 

IV. DERIVATION OF THE STANDARD FORMULAS 

We begin by randomizing the signal parameter, using a general even-order 
chi-squared pdf for a. In particular, we postulate 


f o' a > - I %<!> 


(4-1) 


where g (x) is the pdf defined by Eq. (3-5). 

il 

then 


The characteristic function is 


, /.\ f 1^3 ,a N da 
*o (X) = J e &,(-) 

0 


^i'b' b 


= J e lbXx g,, (x)dx 
0 11 


- a - ibxr 1 , 











a special case of (3-4), namely <^(bA,0). With this model, 

Eu = Ea + N = Mb + N; 

and if tl = kN, the result is the same as assigning to each component signal 
parameter, a , the pdf 


f( V - £ 


y 


with each component independent of the others. In this case, Ea = kb. The 

n 

Swerling cases 1, 2, 3 and 4 correspond to the M-values 1, N, 2 and 2N, 
respectively. 

We note that 

•o tl(1 ~ I 53 = [ (l+bH-b 1 

and abstitute in Eq. (3-18). The resulting expression for the detection 
probability is 


-y 1 


2iri 


lt|=l-< 


,yt 


M-N 


dt 


f(l+b)t-b] h (l-t) 


V 


4-2) 


A basic difference appears now, depending on whether 11 w* or !1<N. We suppose, 
first, that it>N, so the only pole in (4-2) is at t=b/(l+b). We shrink the 
contour to a small circle enclosing this pole and then make the change of 
variable 
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..... . .. 


(4-5) 


the desired expression. 

The. false alarm probability is, of course- 


P FA - V*> 


as it must be for any fixed-threshold model. Note that our pdf, f Q {a), 
approaches 6(a), as h*0, and ^ Q (^) + 1* One can also put 


b = a /It 
o 


and let M+«. In this case f (a) aporoaches 6(a-a ), and 6 (X) exp(ia X). 

o' o o o 

while Pp is transformed into a standard, infinite-series expression for 
P N (y,a o ) (this is Fehlner's series, which is discussed in Section 6). 

For fi=N we have the familiar Swerling case 2 result 


P D “ S N^l+b^ > 


so favored by radar analysts. When H=2!i, we obtain the Swerling case 4 
expression 


p n = jj I Ab 1 C*-2») 

D (1-h>) n £=o 1 1+0 


Whesi M<N we write (4-2) as 


P = e " y -i 
D 2x 


/ ^ 


I ti=l-e ul+b)t-b) t (1-t) 


fa-fi'i 
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and proceed differently. The integrand now has two poles, and we write 


p = P + P 
D r l 2 


(4-7) 


whece the integrands of P^ and P^ are the same as P Q , but the contour for P^ 
is a small circle enclosing the pole at t=b/(l+b), while for it is the 


circle |t|=e. 


The evaluation of P^ is the same as the case just analyzed, which used 


transformation (4-3), up to the point where 


P 1 


_ _L_ 

l+b 1 


ds 


ttt s ■ • -(N-ti) 

/ l+b , s+b s. 

2it:. - | 6 4+b J It.. N 

|s|=e s (1-s) 


This time the expansion is an infinite series: 


-(N-it) , «> 

/■ s+b x /- l+b ,, r ,tHi-l+JL , . .1 

( l+b ) ( "TT> L ( # )(-s/b) , 


fc=0 


which converges since |s|=t' and e can be chosen smaller ‘.nan b. When this 
series is substituted and the integration performed termwise, only the terms 
with £<J1 contribute, since beyond that there are no poles. Theiefore, 


n-m M-l 

p - A+-V V f N-M-l+£ u _ , Jl . y . 

1 ' z l 0 * K 7 Wl?^ ’ 


(4-8) 


where (3-14) has again been employed. 


i . 
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can be handled in a similar way, without changing the variable of 
integration from that used in (4-6) but expanding 


KUMt-bl"" - (- i)" i («-l + *)(«t)‘. 
J £=0 


Again, since |t|=e , convergence can be assured and we get 


P - (-l/b)" ? / «>' _«£_ 

2 \ t-ruJ 2iri . < N-m-j 

£= 0 * D It|=e t X (l-t 


This time the sum stops at £=N-il-l, and the integrals are evaluated using 
(3-14). Altogether, after reordering the suras to get an increasing sequence 
of S-functions, 


M 1 4-k C 

P D ■ P 1 + p 2 - ( - 1/b> <TT> 


f y ( b ' 1 - |1 k- 
I A !H1 - 1 


b >* 


N—11 „ , B , £ 1 


(4-9) 


The usual application of these formulas will be to Swerling cases 1 and 
3, where 11=1 and 11=2 respectively. For these or any small values of 11, 
another expression is more convenient. It is obtained from (4-6) by writing 











(l+b)t-b 



and put 































the complete expression is 


w y > 


(l-x) M 


01 - 1 )! 


(*—) 
v dx' 


11-1 


H 

(d- 


(x-l)y 


x)x 


N-M 


[1 " S N-M (xy)1 


(4-i: 


where 


x 


b 

1+b 


The identity 


df 6X7 


y e xy [1 - S^Cxy)] , 


which is easily verified, makes the evaluation of the derivatives in (4-11) 
relatively straight-forward. 

For Swerling cases 1 and 3 we have: 


N-l-4— 

P D (tl-l) = S^_ 1 (y) + (i±t) e 1+b 


[1 S N-ld+b^ ’ 


(4-1? 


and 


P d (M-2) 


W>-> + (i r )N " 2 ■' lib |a - [i - 


+ i% t 1 -^ 


K><Sb»j 


(4-13 
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By writing 


S (JS-0 
N-ri+b' 


Vife> - e 


by 

l+b 


N-l 


(N-l)! 


(ZL -) 
4+b' 


and substituting into (4-12), it follows that this equation can also be 
written 


yy> 


<¥> 


N-l 


y 

l+b 


[i - S, 


4+b' J 


(4-12) 


In this form, (4-12) is also valid for N*l, in which situation Swerling cases 
1 and 2 coincide. By a similar manipulation, (4-13) can be written 


S N-l (y) + e 7 (l+b)(N--2)! 


N-l 


+ (~) 


N-2 


l+b 


(1-^ 


lfe> 


x t 1 ' S N-l ( I+b )1 (4 ' 13) 

which happens also to be valid for N=2, where Swerling cases 2 and 3 coincide 

We turn now to the dual problem of fluctuating threshold, as in linear 
CFAR detection. The signal parameter, a, is fixed, and the threshold is 
written 

y = cx, 
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where the pdf of x is taken to be g^(x). Thus x corresponds to the sum of L 
random variables, each with the exponential pdf 


g 1 (x) = e x . 


In other words, x is the result of Incoherent summation of L complex samples 
of noise alone, the noise samples having the same variance as the noise 
components of the samples, z^, which contribute to u. We can call these 
noise samples x^, and then 


L x 1 + --* +x L 


y = c l x = cL 


The average of the x is the noise level estimate in CFAR problems, and the 

X 

corresponding threshold multiplier is cL. With this model, the pdf of the 
threshold variable, y, is 


f o (y> * I 8 l { c ) • 


(4-14) 


and the required detection probability is 


P D = / V y,a)f o (y) dy * 

0 


(4-15) 


According to Eq. (3-11), this is obtainable from the solution of the dual 
problem: 
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p D " 1 - I ^(a.y) f D (y)dy 
0 


- 1 - J F 1 -H te,7> g L ( c ) ^ 
0 


This, of course, is the problem we have just solved and we have only to 
make the changes of variable 


N+l-N 

M+L 

y+a 

b-*-c 

in, for example, Eq. (4-2). We note immediately that the case previously 
called '1KN" cannot occur, since the old exponent M-N becomes N+L-l in the 
dual problem. Thus (4-4) provides the desired solution: 


1 - 


(1+c) 


N+L-l 


N+L-l 

l 

£» 0 


(^"L-l \ Z g f a % 

^ Z ; Jt+l-N 4+c y * 


Since ^^(x) vanishes unless m>0, only the terms in this sum having £>K 
survive, and our result can be written in the form 


1 - 


U+c) 


N+L-l 


L-l 

l 

Z=0 


.N+L-l. Z _ , a 

( Ml > c S t+1 ( H? ) 


(4-16) 
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It may be verified that this result is also obtained by a direct 
calculation, starting from (4-15) and using representation (3-8) for P . This 
calculation closely parallels the derivations of Eqs (3-18) and (4-2). 

The false alarm probablity is obtained by putting a=0 in (4-16) and 
noting that S (0)=1: 


P 

"FA 


>J4-T -1 

u+cr 


N-l 


o=0 


,m-L-i. z 
( £ > C » 


(4-17) 


and (4-16) can also be written In the form 


v 

‘D 



c 


(1+c) 


N+L-l 


V w-l-i, _ s 

0 ’ n i N+Z ' - 11 b 

1=0 


£+l'‘l+c' J 


(4-18) 


The probability cf detection when both threshold and signal parameter are 
random variables can, of course, be approached from two directions: we can 
start with random a and fixed y, then averaging over the assumed pdf of y, or 
proceed in the reverse order. It turns out that both methods yield useful 
formulas, not easily transformed into one other. We start with (4-16), 
averaging over the assumed pdf of a: 


f o (a) b 8 mV 


(4-19) 


To carry this out we need the integral 









WIIII|li|llllllllll!ll|iilllllllMI 

. ; 


/ S £+1 (Ax) & a (x) dx 


£ .m ” . 

I h f 6 ^ %< X > dx 

m=0 Q 


f x m 1 r -(1+A)x M-l+m , 
l mf OPTTT 1 e X dx 

m=° o 


-V l • 

(1+Aj* m=0 ® 1 X 


(4-20) 


This expression recurs in several of our formulas, and it should be noted that 
index £ occurs only as the upper limit to the sum. This fact will facilitate 
the recursive calculations of double sums, like (4-21) below. 


When (4-16) is averaged over (4-19), we make the variable change, 
a=bx, to obtain 


P D ~ 1 (1+c) m-L-l ( N+£ 5 c f S £+l ( l+c ) dx ’ 


and now (4-20) provides the desired result 


„ , c" , 1-te >" V ,S!->L- 1 , i V b ,, 

p o • 1 “ -T iffi xi c T+b^> s < mi )c u * < 4 ~ 21 > 

(1+c) £=0 m=0 













Equation (4-21) is valid for all values of It, and it is particularly 
convenient for Swerling cases 1 and 3, where M=1 and 2, respectively, Th 
falsa alarm probability is still given by (4-17). 

Eor 11=1, the ra-sum Js simply 


n-0 


1+b+c 


m 

) = 


1+b+c 

1+c 


£+1. 


1 - ( 


1+b+c 


) 


and (4-21) becomes 


P D = 1 - 


(1+c ) 


N+L-l 


L-l xu , , n ( . £+1 » 

r- N+L-l. £ t , , b . f 

*; 0 ' l [ 


1 - 


(He) 


N+L-l 


N+L-l 

) 

i=K 


f N+L-l. 


l * (--: 


i+b+c 


£+l-N 

) 



The £-sum is now rewritten as the difference between a full sum (£ runnir 
from zero tD N+L-l) and a partial sum, where £ runs from zero to N-l. Th 
final result is (!1=1): 


,l+b N 


N-l 


N-l 


V f W--l, £ ji ,1+bj 


N-l- 


(. 


1+b 


1 
















This formula is convenient because it has N terms, rather than L, and N is 
usually the smaller number. It can also be obtained from (4-12) by averaging 
over y. When N=l, (4—22) becomes simply 


, 1+b 
4+b+c 



(4-23) 


In Swerling case 3 01=2), the m-sum in (4-21) is 

£ 

£ (nri-Du 31 = (1-u) - (£+2)u* + * + (£+l)u^ + ^ l , 

m=0 1 ’ 

where 

_ b 

U ~ l+b-rc 

A calculation which parallels the derivation of (4-22) yields the result 
(M=2): 



(4-24) 


When N=2, this formula reduces to 










u + a+1) 


(4-25) 


P 

~D 


= (■ 


1+b+c 


) 


When!*>N, another formula is more convenient than (4-21). It is obtained 
by starting with (4-4) and randomizing y, thus reversing the order of 
randomization that led to (4-21). The pdf of y will be 





and we get 


P„ = 


(1+b) 1 


II II 

l 

£=0 


V „ ) 


$ Nf£ ^1+b^ g L ( c ) d( c ) 


We change variables (y=cx) and use (4-20) tc find 


P D = ( ‘ 


l+'D 


1+b+c 


) 


(1+b) 1 


M-N 

l 

£=0 


<;) b 


£ 

3=0 


(4 - 26) 


The false alarm probability is still given by (4-17), and the fact that 
(4—26) reduces to (4-17) when b=0 rests upon the useful binomial Identify 


l 


X-l-frnu m 

( _ ) y 


- (l-y ) 3 



(4-27) 


















which is proved in Appendix 1 


When M=N (Swerling 2), (4-26) becomes 


tJ . i- K-l , a 

P « C~~) V ( L “^ )(_§_) 

D U+b+c' v m M l+b+c'' 

m=0 


(4-28) 


For N=I, Swerling cases 1 and 2 coincide, and (4-28) becomes (4-23) in this 
case. Likewise, when N=2, cases 2 and 3 coincide, and (4-28) then reduces 
(4-25), as it should. 

For Swerling case 4, M=2N, (4-26) becomes 


P = (■ 


_l+b_, 

i+b+c J 


(1+b) 


N 


l Q b l < L ><T& 

£=0 * m=0 n 1+b4c 


) 


(4-2S) 


V. GENERALIZATION OF THE FLUCTUATION MODELS 

In Section 4, when a parameter, say y, was randomized, a pdf of the type 



was used, introducing one new parameter. In this Section we generalize by 
using, instead, s "signal-plus-noise" type of pdf, such as 


This generalizes the chi-squared pdf to a non-central chi-squared, and adds 














another parameter. Tne model is equivalent to the postulate y = cx, where x 
has the pdf f^(x,d) as defined by Eq. 3-2. 

We begin by randomizing the signal parameter, a. We put a=bx, where x 
has pdf f M (x,d). In other words 


f 

o 


b f M ( b ’ d ' ) 


(5-1) 


and the characteristic function is 


(X) = E e 
o 


iAbx 


d 

1-ibX 

■d e_ 

M ’ 

(1-ibX)" 


according to (3-4). To apply (3-18), we rpquire 


M 


( M 1 ^ 1 t^ “ 6 ^(l+b)t-b^ exp ^ » 


' (l+b)t-b 


and then 


P D " e 


-y-d 1 


/ e^*" exp [ 


2irx 


td 


(l+b)t-b 


M“N 

t at 


,M„ 


[(l+b)t-b]“(l-t) . (5-2 


It-Ui 


















We confine ourselves to the case M>N and, following the previous analysis, 
make the change of variable (4-3). The new feature is the exponential factor, 
which becomes 


exp [ 


td , ,d s+b. / d . ,1 bd N 

n+bOFb 1 = PXp( s l+b> = exp( T+b ) ex ^e l+b* • 


With the definitions 


d - 


bd 

1+b ’ 


y 


y 

l+b 


we have 


-y-d 


1 

2iri 


e W 


ds 


M,.. V 

s (1-s) 


il-N 

As before, we expand the factor (s+b) , and, this time, evaluate the 

integrals using (3-7). The result is 


(l+b)' 


M-N 


M-N 

l 

&=0 


(**~^) b^ p ( y 

w-ri+b ’ 


bd . 
l+b 


(5-3) 


a direct generalization of (4-4), to which it reduces when d-K). When b=0, 
only the term £=0 remains, and the false alarm probability reduces to the 
expected value 
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The analog of Swerling's case 2 (M*N) is the simple expression 


P 


N'l+b » 1+b' * 


(5-5) 


A result of this kind is to be expected (when M=N), since each complex sample 

z , is now being modeled in the form 
n 


z 

n 


s + w 
n n 


» 


where 


s = Jh (s' + w' ) . 
n n n 


The new "signal components," s', satisfy 









As a consequence, 


z = /l+b (s " + w ") 
n n n 


where 


Elw "I 2 = 2o 2 
n 


and 


? 

3 



Equation (5-5) then follows, since the new effective threshold is y/(l+b). 

In the dual problem, the signal parameter, a, is fixed, while the 
threshold is randomized by the statement: y = cx, when the pdf of x is 
f (x,d). Then 

Li 

■ § £ L ( c > d> <5 ' 6) 

and the detection probability is obtained from (3-11) and (5-3) by making the 
changes of variable 

N+l-N 
M+L 
b+c 
y+a . 
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Parameter d is unchanged, since it plays the same role in both problems. We 
obtain 


(1+c) 


1_ 

N+L-l 


N+L-l m, t « 

l r + ) he 1 p 


/ a x 

£+l-N ( 'l+c * 1+c' s 


(5-7) 


which is a direct generalization of (4-16). Expression (5-7) reduces to 
(4-16) when d=0, since the terms £-0 through £=N-1 will vanish in that limit. 

Although (5-3) contained only P-functions with positive subscripts, (5-7) 
involves negative subscripts as well. The latter are evaluated by means of 
(3-9), but first it is useful to apply (3-9) to all the terms of (5-7), to 
obtain the formula 


(1+c) 


1_ 

N+L-l 


N+L-l . , 

i ( N T > c* ?, 


,_cd _a_. 
N-ri+c » 1+c' 


(5-8) 


Formula (5-8) is the basis of later generalizations, since it is more compact 
than (5-7) and also contains the parameters a and d, in a "natural order.” 

For computation, we split the terms with positive and non-positive 
subscripts, using (3-9) on the latter: 


Pb ' u + =) N+L - 1 


N-l 

V /N+L-l, £ . 

L < £ > c P n-£ (-Si —) 

£=0 X 4+c f 1+c' 


N+L-l „,_ . 

I ( : ) a - ? 


_a_ cd_.J 
£+l-N v l+c ’ l+c' J j 


(5-9) 
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Finally, the first sum in (5-9) is reversed, and the index redefined in the 
second sum to obtain the result 


P 


D 


N 


(1+c) 


N+L-l 


N-l vTLt 1 , A+l „J 

y ( L_ 1 )(i) p (— 

A V L+£ 'V £+l V l+c 
£=0 



+ 


L-l 


l 

£=0 


N+L-l 
• N+£ ' 


£ 

c 


11 ' P £fl ( l+C 



(5-10) 


To obtain the false alarm probability, we put a=0 in (5-10) and note that 

Vi (0 - ’ l - and 


p (SL. 

£+l4+c ’ 


0) “ S £+l<lfe ) 


therefore 


N 


FA 


(1+c) 


N+L-l 


N-l 

I 

£=0 


N+L-l. f l_. 
L+£ } V 


£+1 


S ( c ^ ) 
£+l4+c ; 


(5-11) 


which reduces to (4-17) when d=0. 
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As mentioned in Section 1, the present model describes the performance of 
a linear CFAR system with a non-fluctuating signal and with other signal 
components in the terms used to form the detection threshold. These signal 
components are usually unwanted, and they can strongly modif* the performance 
of a detection system. One obvious effect of these signal contributions is to 
reduce the false alarm probability (by increasing the threshold), which can be 
seen from (5-11). All the terms in that sum are positive, and [cd/(l+c)] 

is a decreasing function of d, since S^(y) is the false alarm probability of a 
system with n samples and threshold y, which decreases with increasing y. The 
linear "CFAR" loses its CFAR property with this change in the statistical 
properties of the noise. It is evident from (5-10) that P^ is reduced by 
these signal components also. 

Since is the probability that u>cx, we have 


1-Pj, = Prob < c] . 

Thus 1-Pjj is the cumulative probability distribution of the "non-central F" 
random variable, p = u/x. Both numerator and denominator of this ratio are 
non-central chi-squared variables, and (5-10) provides a basis for computation 
which is more convenient than the standard series, as given by Springer.(10) 
When d=0, the variable x becomes simply chi-squared, and the resulting 
cumulative probability distribution is a finite sum of elementary functions, 
given by (4-16). It should be noted that all the chi-squared variables 
discussed here are of even order, and this accounts for the relative 
simplicity of the formulas. 

We now extend our results by randomizing a and then d, finally treating 
the case where both a and d are random variables. We begin with a, and 
postulate for this parameter the pdf 
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f 0 (a) 


£ %<!> 


> 


as before. Both sides of (5-8) are multiplied by f Q (a), and integration over 
a is carried out after the change of variable: a=(l+c) x. We get 


P D = 


(1+c) 


n+l-i 


N+L-l 

I 

£=0 


.N+L-l 
^ £ ; 


/ 

0 


> (£i_ X ) o (Us. x ) i+s. dx . 

n-£ 4+c ’ x; b xj b dx 

(5-12) 


The integral is the probability of detection for a fixed threshold and 
fluctuating signal parameter, and hence it is expressible by means of the 
results of Section 4 with the following identification of variables: 


N N-£ 


y 

b 


cd 

* 1+C 


(5-13) 


However, the sum in (5-12) may have to be separated into two parts, where 
N-£5ti in one part, and N-£<31 in the other, unless 11 >N, in which case all terms 
are of the same type. We carry out this procedure only for the case M=1 
(Swerling case 1), before making the assumption that II >N. 

When M=l, we write (5-12) in the form 


P 


D 


( 1+0 


1 _ 

N+L-l 



. N+L-l. 
^ ' 


c Z X(£) 
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(5-14) 


+ 


N+L-l 

l 


£=n 


N+L-l. i 
o ) c 



> 


and use (4-12) to evaluate X(£), and (4-4) for Y(£). First, we recall that 
(4-12) itself can be written 


P D = 


s„(y) + 

N D 


N-l - 


y 

l+b 


tl " 


sO 1 ’ 


(4-12) 


which is valid for all N>1. When the variable substitutions (5-13) are made, 
we find from this version of (4-12) that 


XU) = S„ „( 


cd 


K-t'l+c 


-> + 


N- £-1 


cd 


1+b+c 


[1 ' 


S ( k 
N-Tl+c 


cd 


1+b+c 


)] 


In the corresponding sum over £, we write 


N-l 


I 

£=0 


, N+L-l 


C* X(£) = 


T = s - 1 ^ xcN-i-o 

1=0 x 



to obtain an increasing sequence of S-functions. Then, 
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(1+c) 


r / N+L-l. , 1* £ / Cd v 

\m-i f n 1 h+z n c } ) i+ n+c- 


,1+b+c x £ „ 1+b+c fT _ c , b _cd S1 

+ (___) e [1 ^ +1 (— 


N+L-l ..., . 

+-kr l ( z > c Y(£) • 

(1+C) N 1 £=N 


To obtain Y(£), we use (4-4) with the substitutions (5-13) and '1=1: 


.. 1 " N +* 1-N+£ . V rn , 

Y(£) = ( 1+b+c ) I ( m )( l+£ ) S N+m-£ ( l+b+C ) * 

m-u 


In (5-14), Y(£) enters only in terras where £>N, and also, the S-function is 
zero unless m>£-N. This occurs in only one terra of the ra-sum (the upper 
limit) hence 


!-N+£ 

Y(£) = ^1+b+c 5 S l ( 1+b+c^ 


Of co nT- se 


= e 


and we obtain 

















N-l j N-l 

(l+c )^- 1 I Jo 


,N+L-1.,1*£ /Cd_\ 

K L+£ 'V b £+1'1+c ' 


cd 

1+b+c 


£ f 0 ( L+£ H be ) U S £+ri+c 1+b+c ' 


~ 1+b+c ^ ,N+L-1 W be \ £+ H 

+ e Jo ( »« > ( W {• 


(5-16) 


The terms not involving S-functions in (5-16) can be combined, as follows 


N r 2 ,N+L-1 W 1+b+c/, L r 1 ,N+L-1 W be / +1 
L ( L+£ H ~bT ~ ) + i ( M+9 


N+£ /v 1+b+c ; 


N+ij 1 , T 1 it i £ L L—1 , , L—£ 

V /N+L-l N , 1+b+c N , v /N+L-l N , be \ 

l ^ £ n be ' l ( 9 


£ 1+b+c 


. L .... N+L-l 

<I+sfe> (1+ ^ 


. t(l+b)(l+c)l W ' L 1 
(bc) N_1 (l+b+c) L 
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Therefore finally, 


P = ( i±k) 

D V b ' 


N-l 


L - 

, 1+b N 1+b+c 

( T+b+^ 6 


N-l 


Cl+c) 


N+L-l 


N-l 

I < 

£« 0 


N+L-l 

L+£ 


) e 


£+l'‘l+c ; 


cd 

1+b+c 

e 


(^) £ 


e / b cd '\ 
£+l4+c 1+b+c ; I 


(5-17) 


The corresponding false alarm probability is still given by (5-11), and it can 
be shown directly that (5-17) converges to (5-11) as b+0, and to (4-22) as 
d+0. 

We forego a similar analysis of the Swerliug case 3 problem and return to 
(5-12) with the assumption that M>N. Then (4-4) can be used for all the 
terras, and we obtain 


(1+c) 


N+L-l 


N+L-l 

l 

£=0 


,N+L-l 
1 ; 


( H c - 

±+b+c 


M-N+£ U-N+£ 

) l 

m=0 




ra 


m 


cd 


N+ra-£ 1+b+c 


-) 


(5-18) 
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It is useful to rearrange the sums in (5-18) in order to get an increasing 
sequence of S-functions (to permit their recursive computation), hence we 
write N+m-£ = n+1, and eliminate m. The lower limit on the new index, n, is 
set by the S-functions themselves at n=0, and the upper limit is seen to be 
M-l. In reordering these sums, we note that the original sums over m and l 
are automatically limited by their binomial coefficients, since (”) is zero if 
either k or n-k is negative. Thus, can be written 


P D = 


M-l 

l I 

n=0 < £ 


,N+L-1., M-N+£ . l\ . b ,n 
( £ )( n+l+£-N } X f 


( 


cd 


n+1 v 1+b+c 


) 


(5-19) 


where 


F = 
o 


(1+c) 


M-N-L 


b^d+b+c)™ 




and 

- he 

x = 1+b+c ’ 


The lower limit of the £-sum in (5-19) depends upon n. 
Next, we use the identity 



C-B 

x 


(1+x) 


A-C 



A+£ 

A+B-C 


) 



(5-20) 
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p - c^CL+b)^ M y N ( M-N be * 

D " a+b^)^- 1 »-o * 1+b *= 


C ,N+L-m w l+b “ . cd v 

x i„ ( L+n )( ~> b n+l T+b+c ’ 
n=*0 


(5-21) 


In Swerling case 2 (M-N), the £-sum In (5-21) reduces a single term; and 


, c (1+b) 

D (l+b+c) N+L ~ 1 


N 1 ,N+L-l w l+b n 


V c I ' mj ' J -v±I£) s ( Cd ) 

L ' T +r. « ' VruK+b' * 

n=0 


(5-22) 


The false alarm probability for any M is still given by (5-11), to which 
(5-21) reduces when b=0. 

For Swerling case 4, (5-21) must be used with M=2N. It should be noted 
that (5-21) can be obtained somewhat more directly by averaging (4-4) over the 
pdf, 

f o (y) • I £ l<c - d > - < 5 - 23 > 


for the threshold y. The required ..Legrals, 


/ S f+1 (Xx)f L (x,d) dx, 


are evaluated by writing 
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f L (x,d) - dx P L (x,d), 


then integrating by parts, and using (4-16). 

So far, we nave been discussing extensions of (5-8) in which the signal 
parameter, a, has been randomized. In the dual problem, the parameter d is 
randomized instead. This represents linear CFAR with a fixed signal 
parameter, a, but with fluctuating signal components in the randomized 
threshold. For d we assume the pdf 


f o (d) 



(5-24) 


and average (5-8) with respect to it. The result is 


(1+c) 


N+L-l 


N+L-l hj, T . 

I (“t 1 ) 

2-0 1 


lr D , a . ,1+c 

c ' P N-t* y »l+c* g K ( " 


ch 


y) 


1+c 

ch 


dy (5-25) 


after the change of variable d+(l+c)y/c. We can use (4-16) to evaluate the 
integral, which corresponds to ordinary linear CFAR and fixed signal 
parameter, with the variable changes: 

N+N-A 

L+K 
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a + 


a 

1+c 


c + 


ch 

1+c 


Therefore 


( r\ / a > t 1+C v 1+C j 
J P N-£ (y ’ 1+?' g K (- dh y) “ch dy 


N+K-£-l N+K-£-l . . , m 

_ / 1+c s r ,N+K-£-l w ch 

-*■ M 1_i_l_ ) 


' 1+c+ch' 


l < v >0 S , 

m=N-£ m 1+C 1 


( — - ) 


m /v l+c / m+-£+l-N 1+c+ch 


We assume that K>L to assure the validity of this formula for all £(in 
particular, £=N+L-1). When (4-16) was derived, the parameter there called N 
was non-negative, and we are now applying to a more general case. A review of 
the derivation of (4-16) will show that it is valid as long as N+L-1>0, which 
explains the requirement K>L. This excludes Swerling case 1 type models for 
the fluctuating signal component in the threshold. 

Substituting in (5-25) we get 


/-i , _ L N+L—1 - T T , . , . , £ 

r. _ i (1+c) r ,N+L-l Nr 1+c+ch, 

p °~ 1 'Jo ( * )tc ~ 1 * 


1 ,N+K-£-l u ch . m , a . 

X ( m ^ 1+ch S m+£+l-N ^l+c-;-'^ 
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The sum on m has been allowed to start at m=0, since the S-function will 
vanish in the extra terms. With substitution £-»N+L-l-£, this becomes 


P D " 1 - 


N+L-l K-L N+L-l , 

, 1+c v Y / N+L-l v ,1 1+c , 

^ i £ > L c 1+c+ch J 


' 1+c+ch^ 


K-L+£ . , m 

V ,K-L+4v ,ch» g /_ 

L n K m ' '1+c ' L+m-£ '1+c+ch 

m=0 


) 


(5-26) 


Equation (5-26) can be obtained from (5-]8) by reversing the roles of the 
original random variables, complimenting the probability, and making the 
parameter substitutions 


N+L 

L+N 

ll+K 


b-*h 

d+a 

c+l/c 


(5-27) 


A form analogous to ’5-21) can be obtained from (5-26), either by a parallel 
calculation or by making the variable substitutions (5-29) in (5-21). Either 
way, the result is 


1 - 


K /i.u,\ N + L - K 

c (1+h) 


(1+c+ch) 


N+L-l 


K—L , r . . SL 

i c'; l ki h 


£=0 


1+c+ch 


) 


x 


L-l+.l 


n=0 


.N+L-l+JL , ,.n c 

( s+n X'* 10 s 


< 


n+1 1+c+ch 


) 


(5-28) 
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It can be seen that (5-28) reduces to (4-16) when h=0, since then only the 
tern £=0 survives. 


The false alarm probability is obtained from (5-28) by putting a=0, which 
has only the effect of replacing the S-function by unity. But then. 


L-l+£ NfL-l+£ 

i <^; m )<c-*h)”, < c4ch) -» i ( ‘ ,+L : i+i )( C ^h)" 

n=0 U " n-N 


(1+c+ch) 

(c+ch) f 


N+L-l+£ 


N-l 

1 r ,N+L-1+£ W , ,.n 

- s l ( )(c+ch) , 

(c+ch) n=0 


(5-29) 


I ( K : L )h"= (l+h) K ' L . 
£='0 * 


The result is that the first term of (5-29) eventually cancels the unity in 
(5-28) and the false alarm probability simplifies to 


/, ,u\I j —K K—L £ 

p sx —L Lill - y ( K —) x 

FA N+T .-1 4 v 5 ' a 


FA , S N+L-1 on i 1+c+ch 

(1+c+ch) £=»0 


r ,N+L-l+£v, ,^n 

x ) ( )(c+ch) 


(5-30) 
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This expression reduces, of course, to (4-17) when h=*0. 

The final situation, in which both of the signal parameters, a and d, 
are randomized is analyzed first by averaging (5-28) with respect to the usual 
pdf for a: 


£ o (a> = k %<!> 


In this way, a formula valid for all M is obtained, although the constraint 
K>L is still in force. Only the S-function in (5-28) depends on a, and the 
use of (4-20) gives us 


/ S n+l ( l+ c +ch ) %^b^ d( b ) 


, 1+c-tch > M r b 

^1+b+c+ctv m 1+b+c+clr 

m=0 


(5-31) 


When (5-31) is used to average (5-28), the result is 


n N+ l ~K ,, , , , -M-N-L+l 
D _ , c (1+h) (1+c+ch) 

D M 

(1+b+c+ch) 


K ~L v_ T U l L-l+£ 

l J 0 - 
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X 


m=0 


,11-1+nK > b ' 
m '^1+b+c+ch' 


m 


(5-32) 


and the probability of false alarm is still given by (5-30). 

To reverse the constraints, we can obtain a formula valid for all K and 
11 >N by averaging (5-21) with respect to the pdf of d, namely (5-24). The 
resulting detection formula is 


P 


’D 


N-l 

c 


(l + b) N ^Wc ) K ~ N ' L+1 

(1+b+c+ch)^ 


x 


x 


11-N 


l 

£=0 



)( 


be 


1+b+c 


) 


N-l+£ 

I 


n=0 




X 


X 


n 


I 

m=0 


•K-l+nu, ch 
m ^ 1+b+c+ch 


m 

) 


(5-33) 


Expression (5-33) reduces to (4-26) when h=0, as does (5-21) in the 
corresponding limit, d=0. It is possible to show that (5-32) and (5-33) are 
equivalent when both K>L and 11>N. 

The false alarm probability corresponds to b=0: 
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X 


c N- 1 (i +c) K-N-L + 1 N-l n 

-- -K L K L+n n c J 


(1+c+ch) 


f K-l+m v ch u 
k m n l+c+ch ; 


(5-34) 


We did not work out the detection and false alarm probabilities for the 
analogous fixed-signal case, with K<L, but (5-34) is identical to the result 
of averaging (5-11) over the pdf (5-24) for d. 

Formulas (5-32) and (5-33) generalize some results of Rickard and 
Dillard/ 11 ^ and also Wishner’s analysisof the normalized periodogram 
detector. 

VI * RELATION TO OTHER FORMS 

A striking feature of the finite sum detection formulas is the diversity 
of equivalent forms they can assume. Solutions to the same problem, such as 
(4-21) and (4-26), deduced from different starting points can be reconciled 
only with difficulty. In this section, we relate a few of our results to 
quite different formulas which have appeared in the literature. 

Consider Eq. (4-16) which describes linear CFAR with a fixed signal 
parameter, a. It is 


1 - P =--- 

D a+c) S+L-) 


which can be written 


,N+L-1, £ 

N+£ )C S £+l ^1+c 5 » 


( 6 - 1 ) 
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N 

•) 


(1 


V 


F N a-D 


L-l 

l T (L-l, £) 
£=0 


( 6 - 2 ) 


where 


T n (L,£) 


1 ,N+L. 2 C 

,L C S J1 + I1+C ) 

(l+c; 


(6-3) 


How f n (L) can be rewritten as a sum of terms, each of which depends only on £, 
and not also on L, as do the terms T^(L, £). To accomplish this we write 

F »‘ L > ■ X t n ( l >« ■ V L - U + V L -« 

£=0 £=0 


L-l 

- T (L,L) + l [T (L,£) - T (L-l,£)] + F„(L-1), 
N £=0 N N N 


or 


F n (L) = A N (L) + F N (L-1), 


(6-4) 



where 


L-l 


AfL) = T„(L,L) + l [T„(L,£) - T„(L-1,£)3 
h N £=0 ” 
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(6-5) 



ill 











for L>1. Iterating (6-4) we obtain 


t N (L) = A N (L) + a n ( L_ 1} + *** +A N ( 1> + V 0) * 


Defining 


V 0) = F n ( 0) = I Kt (0,0) , 


( 6 - 6 ) 


we can write 


V L > 



(6-7) 


which has the desired characteristic. 
As applied to (6—2), we have 




1 

(l+e) L 


,N+L. * c , N 


( 6 - 8 ) 


where 



One finds that 
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-u r ,N4L-1. (cu) * 
i=o 


(6-9) 


which holds alsc for L=0. Tnerefore 


F„(L) = e U l 


r (cu) m 

i V fi+T--1 1 ml 


£=0 (14c ) 1 ra=0 llHr ' 1 m! 


and L now appears only as the upper limit of the outer sum. 


me finite nt-sum is a special case of the confluent hypergeonetric 

/12\ 

function; in fact it is a Laguerre polynomial. The basic definition is 


L (k) (x) = * ( k4n } 

m lx) i 'k+a^ n! 


( 6 - 10 ) 


hence 


c n -v _ ~ u r 1 - (N-l), , 

f n (L) - e l -- £ L £ (-cu) , 


£=0 (14c) 


and finally 


„ K "iXT L ~ l 

**-!-&> * ^ l 


£=0 (14c) 


,(N-1) , ca . 
t 14c } 


( 6 - 11 ) 






The Laguerre polynomials can be computed from the second-order recursion 

. ( 12 ) 

l elation 


L< k) (-x) = (2 + L ( £ k >(-x) ~ (1 + L ( £ k > (-x) 


( 6 - 12 ) 


together with 

L^ k) (-x) = 1 

(6-13) 

Lj k) (-x) » k+l+x . 


(13) 

The particular case, N=l, of formula (6-11) was given by Finn, whose 
L £ (x) equals it L £ ^(-x), in our notation. 

It is interesting that the probability of false alarm obtained from 
(6-11) appears in the fora 


P = 1 - 
FA v l+c' 


L-l 

r 

L 

£= 0 


<T> 


l 

£ * 

(1+c) 


(6-14) 


which can be transformed into (4-17) by the use of (4-27). 

Formula (4-4) can be treated the same way, when M>'J. According to (4-4) 


Pn = K cm) 

































hence 


bu 


* ,„x uu -u (K-l)! , (N) , N 

V K) (l+b) K 6 (N+K-l)! L K-1 ( " bu) 


Finally, 


F N (K) - s (U) + I V K > 

k-l 


„ , „ , , N -u £ (k-l)! 1 t (N) , , s 

= S N (u) + bu e L k -1 ( ' bu) > 


M_ y _1 k! 1 (N) by_ 

k=0 (N+k)! (l+b) k k l+b 

(6-16) 

This time, the false alarm probability (b=0) is simply 

P FA ' V y) > 

the same as (4-5). 

We demonstrate another kind of equivalence now, reducing an infinite 
series solution for the fixed-signal, linear CFAR problem to our finite sum 


and therefore 


„ s (-Z-) + — (-Z-) 

Vl+b ; l+b 'l+b'' 


_ _I_ 

l+b 
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(4-16). The former solution is obtained by averaging Fehlners' series: 


P N (y»a) 





(6-17) 


over the pdf 


£ „<y> 



fcr the threshold. Series (6-17) can be obtained from (3-7) by expanding the 
factor exp(a/t) in a power series and integrating term by term. These 
integrals are performed by using (3-14). 

According to (4-20), 
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The last step in the derivation of (6-18) is an application of (4-27). 
Another transformation yields 
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When (6-19) is used to average (6-17), the result is 
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But 
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in terms of the confluent hypergeometric function/ 15 ^ Applying Kummer's 
first transformation/ 16 ^ 
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The new hypergeometric function is simply a polynomial, since its first 
argument, m+l-L, is a non-positive integer. Altogether, 
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Next, we replace is by the new index n = m+4 , and find 
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which is -i&). Similarly by averaging (6-17) ever a, u-ing the usual pdf 
with M>N, (4-4) can be derived. 
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APPENDIX 1 


The binomial identity 
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m=0 
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( „ > x 


(l-x) n 



(A-l) 


is proved by taking the factor (l-x) n inside the sum and expanding again, as 
follows: 
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In the last step, we have rearranged the sum in increasing powers of x. now 
the m-sum is 
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■ 0 j 


(A-2) 


according to a standard binomial identity, which completes the proof. The 
standard identity used here may be written in the form 


l (g)( £ "^ _S )(-l) S = (-lA k ”*> , (A-3) 

s=0 


which is valid when k>£, and is proved by use of the obvious identity 


(l+x) k (l+x)“ £ = (l+x) k 1 . 

Each factor on the left is expanded in a binomial series and the resul 
double sum reordered in powers of x. The desired result follows by equating 
to the expansion of the right side. 

To prove (5-20), we need an intermediate result. Consider the identity 
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[’+x(l+y}] m = (!■*)“ (1 + yg-) 
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and expand the left side and the right-most factor on the right side: 
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Explicit, summation limits are not used because the sums are limited 
automatically by the binomial coefficients, as mentioned in connection with 

g 

Eq. (5-19). Now we multiply both sides by (1+y) and expand some more. Th 
left side becomes 
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and the right side becomes 
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we introduce the new index, k=n+j, and eliminate n, so that the last double 
sum expression becomes 
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Equating powers of y, we obtain the desired result: 
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It is necessary only that m, s and k be non-negative, and the ranges of the 
sums over £ aid j, as limited by the binomial coefficients, will depend on the 
relative magnitudes of m, s and k. 

Now consider (5-20): 
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Using (A-5) of che left side we get 
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and applied to the right: 
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In the last step we replaced j by A-C+j, and this completes the proof since 
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